
clear;
clc
w=0:0.01:2*pi;
j=sqrt(-1);

figure(1);
for p=1:4
H(p,:)=rootlocus_AB(exp(j*w),p);
plot(H(p,:));
hold on;
end
xlabel('Real Axis'),ylabel('Imaginary Axis'),title('rootlocus_{ABF}'),grid on
saveas(1,'root-locus-ABF.png');

figure(2);
for p=3:5
H(p-2,:)=rootlocus_AM(exp(j*w),p);
plot(H(p-2,:));
hold on;
end
xlabel('Real Axis'),ylabel('Imaginary Axis'),title('rootlocus_{AMF}'),grid on
saveas(2,'root-locus-AMF.png');

figure(3);
for p=1:4
H(p,:)=rootlocus_BD(exp(j*w),p);
plot(H(p,:));
hold on;
end
xlabel('Real Axis'),ylabel('Imaginary Axis'),title('rootlocus_{BDF}'),grid on
saveas(3,'root-locus-BDF.png');